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Motivated by the tension between the Higgs mass and muon ^ — 2 in minimal supersymmetric 
standard model (MSSM), we analyze the muon g — 2 in supersymmertic B — L extension of the 
standard model (BLSSM) with inverse seesaw mechanism. In this model, the Higgs mass receives 
extra important radiative corrections proportional to large neutrino Yukawa coupling. We point 
out that muon g — 2 also gets signihcant contribution, due to the constructive interferences of light 
neutralino effects. The light neutralinos are typically the MSSM Bino like and the supersymmetric 
partner of U (1)b-l gauge boson (Hffno). We show that with universal soft supersymmetry breaking 
terms, the muon g — 2 resides within 2cr of the measured value, namely ~ 17 x 10“^°, with Higgs 
mass equal to 125 GeV . 

PACS numbers: 


The Standard Model (SM) prediction for the anoma¬ 
lous magnetic moment of the muon, = (^ — 2)^/2 
(hereafter muon g — 2) has a discrepancy with the exper¬ 
imental results: 

Aa^ = a“P - a™ = (28.7 ± 8) x 10“^° (Icr). (1) 

This discrepancy has survived after performing highly ac¬ 
curate theoretical calculations [1] within the SM frame¬ 
work and experimental analyses [2]; and hence, it 
can be resolved or ameliorated by contributions from 
new physics beyond the SM (BSM). If supersymmetry 
(SUSY), as one of the forefront candidates for the BSM 
physics, is a solution to the muon g — 2, the SUSY par¬ 
ticles, namely, smuon and weak gaugino (Bino or Wino) 
masses should be around a few hundred GeV, in order to 
utilize the supersymmetric contributions [3]. 

However, the observation of the Higgs boson of mass 
about 125 GeV requires rather heavy sparticle spectrum 
within the MSSM framework, and it results in a strong 
tension in simultaneous resolution for both the 125 GeV 
Higgs boson and the muon g — 2 problem since SUSY con¬ 
tributions to muon ^ — 2 is suppressed by the heavy spec¬ 
trum. Non-universality in gaugino and/or scalar masses 
may remove this tension [4], nevertheless in this case 
SUSY models will have plenty of free parameters and 
will lose their productivity. 

In this article we show that this tension can be alle¬ 
viated in the U{1)b-l extended Supersymmetric Stan¬ 
dard Model (BLSSM), which is one of the interesting non- 
minimal realizations of supersymmetry. The BLSSM is 
also well motivated by the established existence of non¬ 
zero neutrino masses [5]. It turns out that, in this class 
of model, the scale of B — L symmetry breaking can be 
related to the scale of SUSY breaking [6], therefore, a 
TeV scale type I or inverse seesaw mechanism can be 
naturally implemented [7]. The one-loop radiative cor¬ 
rections to the SM-like Higgs boson mass, due to the 
right-handed (s)neutrinos in BLSSM, with inverse see¬ 
saw mechanism, provide new contribution to the Higgs 
boson mass in addition to the stop sector of MSSM [8]. 
Thus the lower bound imposed by Higgs mass on the uni¬ 


versal gaugino soft masses mi /2 is reduced, which makes 
possible to find solutions with the light weak gauginos. 
Moreover, in BLSSM with inverse seesaw the g—2 may re¬ 
ceive new contributions, in addition to the usual MSSM 
ones, due to the extension of the neutralino sector by 
SM singlet {B — L) Higgsino and B'-'mo and also due to 
the possibility that one of the right-handed sneutrinos 
is light (due to large mixing between right-handed sneu¬ 
trinos and right-handed anti-sneutrinos). We emphasize 
that within BLSSM with inverse seesaw the g — 2 resides 
within 2<j of the measured value. We show that in this 
model, even if we are restricted to the universal bound¬ 
ary conditions at the GUT scale, we can account for the 
Higgs mass and keep the lightest neutralino to be as light 
as about 200 GeV. 

The BLSSM is based on the gauge group SU{3)c x 
SU{2 )l X U{1)y X U{1)b-l^ One can implement the 
seesaw mechanism by adding the right-handed neutrino 
one per family as required by the anomaly cancellation. 
In order to count for the right-handed neutrino contribu¬ 
tions, we employ the inverse seesaw model in our study. 
The BLSSM with inverse seesaw includes, in addition to 
the MSSM particle contents, two SM singlet chiral Higgs 
superfields xi ^2 and three sets of SM singlet chiral super¬ 
fields Ui, 51-, 52- [9]. The Superpotential of this model is 
given by [10] 

w = -gr]XiX2 PgHuHd Pgshh -YgdqHd 
— Y^elfld -\-YuUqHu + z)xi <§2 + W z)/(2) 

The definition of the parameters appear in W, the corre¬ 
sponding soft SUSY breaking terms and the details of the 
associate spectrum can be found in Refs. [9, 10]. Note 
that the U{1)y and U{1)b-l gauge kinetic mixing can be 
absorbed in the covariant derivative redefinition. There¬ 
fore, this mixing provides another source of coupling be¬ 
tween MSSM and B — L sectors. Here, we will focus only 
on the particles involved in the g — 2 loops, namely light 
neutralino, smuon and chargino, sneutrino. 

Considering the additional SM singlet fields of BLSSM 
with inverse seesaw mentioned above the B — L exten- 
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sion may modify only the neutral sectors of the MSSM 
only, and hence the chargino and slepton mass matri¬ 
ces remain intact. The 7x7 neutralino mass matrix, in 

the basis: 5', xi, X 2 ^, can be found in 

Ref. [10]. One can easily show that depending on the 
ratio of the couplings gi and Qbl^ the lightest neutralino 
could be 5-ino {B) or B'-ino {B') like. It is worth noting 
that in order to account for the Higgs mass, mo should 
be of order TeV, thus /r and which are determined 
from electroweak and B — L symmetry breaking condi¬ 
tions, respectively, are also of order TeV scale. Hence, 
the lightest eigenstates of neutralinos are mostly formed 
by the neutral gauginos, since the neutral Higgsinos are 
quite heavy. Now we turn to the sneutrino mass matrix. 
If we write t)l,r and S 2 as t)l,r = + icFL^R) and 

S 2 = + '^^ 5)5 then we can get the CP-odd/even 

sneutrinos matrices as given in [10]. 

The supersymmetric contributions to in BLSSM 
can be split into neutralino and chargino parts as for 
the MSSM [3], 
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where Ufi are the rotation matrices which diago¬ 

nalize the muon and smuon mass matrices respectively, 
while Ufyi and Uf^n diagonalize the CP-odd and CP-even 
sneutrino mass matrices. Note that one can neglect the 
mixing between slepton families, and consider the smuon 
mass matrix as 2 x 2 matrix separately from the first 
and third families. In addition, the mixing between two 
smuons is proportional to the Yukawa coupling associ¬ 
ated with muon, which is of the order ^ 10“^, and hence 
left- and right-handed smuons are approximately match 
with the mass eigenstates, and hence, the rotation ma¬ 
trix for the smuons, Ufi can be set to unity in a good 
approximation. A similar discussion holds for the muon 
mass matrix diagonalized by 

As seen from Eqs.(5,6), the Bino contribution is in a 
similar form as obtained in the MSSM, but in BLSSM it 
has an enhancement by the gauge mixing between U{1)y 
and U{1)b-l characterized by the coupling g. It is worth 
emphasizing the contribution from H'—ino {B'). It con¬ 
tributes to through interactions with muon governed 
by H — L gauge group. In addition, its contribution ex¬ 
hibits an enhancement with the gauge kinetic mixing. 
Moreover, since it is allowed to be as light as Bino, and 
even lighter, the lightest neutralino can be formed to be 
mostly B' OT B — & mixing. Thus, one can expect its 
contribution to be comparable with that from Bino in 
BLSSM; i.e. Nai ~ numerically. We also present the 
contribution from the Higgsino component of the Neu¬ 
tralino; however, due to the smallness of Y^, its contri¬ 
bution is strongly suppressed. 

Similarly Eqs.(7,8) reveal the contribution from the 
chargino expressed in terms of CP-odd and CP-even 
sneutrino sectors separately. Note that Eqs.(7,8) hold 
approximately, and one can combine these two sectors 
in the case of strong mixing between them by summing 
over j = 1,..., 9. The contribution denoted by are 
mostly suppressed because of Y^. On the other hand, 
is expected to dominate in the chargino contribu¬ 
tion to a^. It arises from the interactions between muon 
and sneutrinos through SU(2) interactions as shown in 
the first and third terms in Eq.(8). The remaining terms 
are from the Yukawa interactions, and they cannot be ne¬ 
glected, since Yy is allowed to be of the order 0(1) when 
one employs the inverse seesaw mechanism in BLSSM. 

In this regard, the parameters relevant to muon ^ — 2 in 
BLSSM with inverse seesaw can be listed as M^,, 

and Yy in addition to those present in MSSM, where 
M^, is the soft SUSY breaking mass term for the H', 
rriyi and are the masses of the lightest CP-even 

and CP-odd sneutrinos respectively. Besides the mass 
eigenstates, g, which quantifies the kinetic gauge mixing 
between U{1)y and U{1)b-l^ and Dirac Yukawa cou¬ 
pling associated with the neutrinos denoted by Yy are 
also effective in supersymmetric contributions to a^. In 
scanning the parameter space, we have employed the 
Metropolis-Hastings algorithm as described in [11], and 
used SARAH [12] and SPheno [13] for the numerical re¬ 
sults. The data points collected all satisfy the require- 
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FIG. 1: Aa^ as function of m^o (top) and ?ti~± (bottom). 


ment of radiative electroweak symmetry breaking. After 
collecting the data, we impose the mass bounds on all the 
particles [14]. We have employed the Higgs mass bound 
as 123 GeV < rrih < 127 GeV [15, 16], where we take 
into account about 2 GeV uncertainty in Higgs boson 
mass due to the theoretical uncertainties in calculation 
of the minimum of the scalar potential, and the experi¬ 
mental uncertainties in rrit and as. We also employ the 
gluino mass bound: rrig > 1 TeV [17] and the neutral 
gauge boson Z' mass bound: Mz' > 2.5 TeV [18]. 

We present our results first in Fig. 1, where we plot the 
muon ^ — 2 as function of lightest neutralino (top) and 
lightest chargino (bottom). The green points satisfy the 
mass bounds. As can be seen from the panel Aa^ — m^o, 
muon g — 2 puts a sharp upper bound on the neutralino 
as m^o < 400 GeV. Similarly, Aa^ — exhibits a 
sharp bound on the lightest chargino mass at about 600 
GeV. Also the chargino cannot be lighter than about 500 
GeV because of the mass bounds we employed in our 
analysis. Therefore, the sneutrino-chargino contribution 
is rather suppressed by the chargino mass. Both planes 
indicate that the best result for muon ^ — 2 is found about 


17 X 10 consistently with the constraints. 



FIG. 2: The correlation between the masses of B' and H(top) 
and Aag versus tan /3 (bottom). 


In Fig. 2 we display the correlation between the masses 
of the light neutralinos; B' and B (top) and also the 
dependence of Aa^ on tan^^ (bottom). The color coding 
is the same as Figure 1 except the plane M^, — Mi in 
which the yellow points are a subset of green and they 
satisfy the muon g — 2 condition within 2cr. The unit line 
in this plane indicates the region with M^, = Mi. As 
can be seen from the M^, — Mi plane, M^, is lighter than 
Ml. Hence, Bino and H'—ino dominate in muon g — 2 
results through the mixing of neutralinos. According to 
the results, M^, is bounded in the interval ^ 200 — 250 
GeV, while Mi in ^ 250 — 360 GeV. One can obtain 
some significant contributions to muon ^ — 2 in such a 
narrow mass intervals for B' and B with the help of tan p 
enhancement as shown in the plot of Aa^ — tan P that 
muon g — 2 within 2cr requires tan /3 > 45. 

It is worth noting that the contributions from M^, and 
M^ is further enhanced by their gauge kinetic mixing in 
smuon-neutralino channel as given in Eqs.(5, 6). In Fig. 
3 we show the dependence of muon ^ — 2 on the gauge 
kinetic mixing characterized with g. The color coding is 
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the same as Figure 1. The figure indicates that ^ ^ — 0.1 
is favorably for enhancing Aa^. 
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FIG. 3: Aa^ as function of g. The color coding is the same 
as Figure 1. 

In conclusion, we have found that the supersymmetric 
contribution to muon ^ — 2 in BLSSM with inverse see¬ 
saw mostly relies on the light Bino with Mi > 250 GeV, 


and even lighter B' of mass about 180 GeV. In this case 
the lightest neutralino mass eigenstate is either mostly B 
or a linear superposition of B' and B. The tan^d depen¬ 
dence of Aa^ is also represented, and we have found that 
tan P needs to be > 45 in order to raise the supersymmet¬ 
ric contribution such that muon g — 2 results satisfy the 
measurement within 2cr. In addition, the kinetic gauge 
mixing coupling g enhance the supersymmetric contribu¬ 
tion to muon g — 2. We have revealed in our analysis 
that muon g — 2 favors g —0.1. The gauge kinetic 
mixing affects the contribution from smuon-neutralino 
channel, and with its enhancement, the smuons can have 
the mass about a TeV, and hence the higher mo values 
can still yield results compatible with muon — 2 as well 
as consistent with the 125 GeV Higgs boson mass. 
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